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Abstract 

Evanescent modes in complete sonic crystals (SC) and SC with point defects are reported both 
theoretically and experimentally in this paper. Plane wave expansion (PWE) and, in general, u(k) 
methods have been used to calculate band structures showing gaps that have been interpreted as 
ranges of frequencies where no real k exists. In this work, we extend PWE to solve the complex 
k(u>) problem applied to SC, introducing the supercell approximation for studying one vacancy. 
Explicit matrix formulation of the equations is given. This k(u) method enables the calculation of 
complex band structures, as well as enabling an analysis of the propagating modes related with real 
values of the function k(u)), and the evanescent modes related with imaginary values of k(u). This 
paper shows theoretical results and experimental evidences of the evanescent behavior of modes 
inside the SC band gap. Experimental data and numerical results using the finite elements method 
are in very good agreement with the predictions obtained using the k{oj) method. 

PACS numbers: 43.20. +g, 43.35. +d, 63.20.D-, 63.20.Pw 
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I. INTRODUCTION 



The propagation of scalar waves inside periodic structures has been receiving growing 
interest in recent years. A great effort has been made to understand the physics of these 
systems since the acoustical properties of a periodic sculpture by Eusebio Sempere were 
measured.- 

Phononic crystals (PCs) consist of an inhomogeneous periodic distribution of elastic ma- 
terials embedded in other elastic materials with different properties.- 1 ^ These systems are 
extensions of the photonic crystals 4,5 used for the propagation of elastic waves through peri- 
odic elastic structures. If one of the elastic materials is a fluid medium, then PCs are called 
sonic crystals (SC). Several studies discuss the similarities and differences between these 
periodic systems.-^- 

The periodicity of these systems is introduced in the solution of the wave equation by 
means of Bloch's theorem. This solution leads to the phenomenon of band gaps (BGs): 
frequency regimes where waves do not propagate through the crystal. Traditionally, wave 
propagation inside such systems was analyzed by means of the band structures. Plane wave 
expansion (PWE) 8 transforms the wave equation into an eigenvalue problem that can be 
solved for each Bloch vector, k, in the irreducible first Brillouin zone; and so obtaining 
the eigenfrequencies u(k) that constitute the band structures. In the case of SCs, it has 
been proven that eigenfrequencies for an arbitrary crystal structure and an arbitrary filling 
fraction" are real values. A great number of applications based on SCs are explained by the 
existence of BGs: acoustic filters;— acoustic barriers;— or wave guides.— 1 ^ 

Propagating waves inside a periodic media represent a set of solutions to the wave equa- 
tion that satisfy the translational symmetry, and these are characterized by the transmission 
bands in the PWE method. However, where the translational symmetry is broken, finite 
periodic media or periodic media with point defects, can support the well known evanescent 
modes characterized by a complex wave number, k.— Recent experimental results^ show 
measurements of the sound levels recorded inside a point defect and behind an SC. These 
authors observed that this level is higher inside the cavity than behind the crystal. This 
fact clearly shows both the generation of a trapping mode (i.e. localized mode) inside the 
point defect and its evanescent behavior outside the vacancy. Some authors in the elec- 
tromagnetic regime have measured the evanescent modes in photonic crystals and revealed 
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multi-exponential decay.— 

Several extensions of the PWE method have been used to analyze the propagation of 
sound through periodic systems in different situations; for example, crystals with point 
defects have been analyzed with PWE using the supercell approximation.—*^ The same 
methodology has been used to analyze the influence of the following: constituent materials, 
plate thickness, and the geometry of the array on the band structure in two dimensional 
(2D) phononic crystal plates.— However, these u(k) methods interpret the BG as frequency 
regimes where no real k exists. Therefore, these methods can only be used to study and 
characterize propagating modes. 

We have been motivated by the work of Hsue et al.,— in which the PWE was extended 
for the case of photonic crystals to calculate the complex k in a 2D isotropic and in general 
3D anisotropic cases. In this paper we show the extended plane wave expansion (EPWE) for 
the case of 2D SCs. The aim is to obtain the band structures using the inverse expression 
k(oj), and with a possibly complex k. Recent works show the calculation of complex band 
structures for phononic crystals.— 1 ^ In the present work we show the explicit matrix for- 
mulation and the approximation of supercell for analyzing the complex relation dispersion 
of SCs. The extension of the methodology enables us to characterize the evanescent and 
propagating modes in complete SCs, as well as in SCs with point defects. 

In this paper we present novel measurements of the pressure in the space between rows 
inside an SC. We have developed a 3D computer-controlled automatic positioning system to- 
gether with an automatized acquisition system, called 3DReAMS (3D Robotized e-Acoustic 
Measurement System). This system enables the pressure field in trajectories inside a crys- 
tal to be measured, and we have consequently analyzed the decay of the evanescent modes 
throughout an SC. The imaginary part of the wave number of the evanescent modes can be 
obtained experimentally with the measurements taken by 3DReAMS. These data represent 
the experimental confirmation of the analytical results obtained by the EPWE, as well as 
an experimental analysis of propagating and evanescent modes in an SC. 

The paper is organized as follows. Section [III summarizes the main ingredients of the 
PWE for 2D SCs with the explicit matrix formulation of the problem. In Section II III we 
extend the PWE to the EPWE to solve the eigenvalue problem k(oj). We show the matrix 
formulation, as well as the EPWE, together with the supercell approximation for studying 
the complex band structures of 2D SC with point defects. In Section [TV] the complex band 
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structures of an SC of PVC cylinders embedded in air are obtained with EPWE for a 2D 
SC with, and without, point defects. Experimental results validating the predictions of the 
EPWE for the evanescent and propagating modes are shown in Section |V] Finally, the work 
is summarized in Section IVII 



II. PLANE WAVE METHOD 



Propagation of sound is described by the equation 



ld2p vfivp) (1) 



pc 2 dt 2 \p 

where c is the sound velocity, p is the density of the medium, and p is the pressure. 

In this paper we consider a system composed of an array of straight, infinite cylinders 
made of an isotropic solid A, embedded in an acoustic isotropic background B. There is 
translational invariance in direction z parallel to the cylinders' axis; and the system has a 
2D periodicity in the transverse plane. By making use of this periodicity, we can expand 
the properties of the medium in the Fourier series, 

°=^ = Y.°i(G)e' S ', (2) 

G 

^ = ^ = E^)^- ( 3 ) 
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G is the 2D reciprocal-lattice vector and B(r) = p(r)c(f) 2 is the bulk modulus. For the 
pressure p we use the Bloch theorem and harmonic temporal dependence, 
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where (3 = (<J,rj), and F(G) is the structure factor. For a circular cross section of radius r, 
the structure factor is 



A uc is the area of the unit cell, A cy \ is the area of the cylinder, and J± is the Bessel function 
of the first kind of order 1. 

Using equations (J2J), (J3]), fll]) and ([T]) we obtain^ 

^2 ((k + G)a k (G - G>){k + G>) - u 2 V (G - G')) P%{G') = 0. (7) 
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For G taking all the possible values, Equation ([7]) constitutes a set of linear, homogeneous 
equations for the eigenvectors Vk(6) anc ^ eigenfrequencies co(k). We obtain the band struc- 
tures when k scans the area of the irreducible region of the first Brillouin zone. 
Equation (J7J) can be expressed by the matrix formulation below 



where i=l,2,3. The matrices Tj, £ and Q are defined as 
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The explicit matrix formulation is shown as follows: 
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where G = (G x , G y , G x ). To solve ([HD we must truncate the matrices. If we chose m = n = 
(-M, . . . , M), the size of the previous matrices is N x N = (2M + 1) x (2M + 1). NxN 
is usually the number of plane waves used in the calculation. 



By solving the system given in (JSJ) for each Bloch vector in the irreducible area of the 
first Brillouin zone, we obtain N x N eigenvalues, u 2 , which can be used to represent the 
band structures, ui{k). 

III. EXTENDED PLANE WAVE METHOD 

In the u(k) formulation, the existence of BG is indicated by the absence of bands in deter- 
mined ranges of frequencies. However, BG could be understood by means of the evanescent 
behavior of the internal modes. This interpretation was predicted by some authors^ when 
approximating the second band near the BG by expanding u(k) to powers of k around the 
edge k = n/a, being a the lattice constant of the array. These authors claimed that as the 
BG is traversed, the exponential decay grows as the frequency nears the center of the BG. At 
a given frequency u inside the BG, the evanescent wave is characterized by a complex value 
of its wave number k(u) and which the imaginary part characterizes as the exponential-like 
decay of the mode. In this section, we extend the previous PWE to the EP WE to obtain 
k(co) and with a possibly imaginary k. 

From Equation (JSJ) we define the next vector, 

= sr,p. (14) 

With this definition we can reformulate the eigenvalue problem (JSJ) as the equation system 

3 

U 2 QP = J2 T i®i- ( 15 ) 
i=l 

To obtain an eigenvalue problem for k(u), we write k = ka, where a is a unit vector. Then 
f fTUj) can be written as 

r, = r° + kaj, (16) 
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where I is the identity matrix, and 
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Equation (JHI) can then be written as 



^ - Eli r°sr° o 
- Eli sr° / 







H 









Eli r?E«i / 

ELlSok o 









I ^ ) 



(19) 



where $' = J2i=i a i^i- 

Equation ( fT9l) represents a generalized eigenvalue problem with 2N eigenvalues k, and 
possibly complex numbers for each frequency. Complex band structures have been calculated 
for the incidence direction characterized by vector a by solving the previous eigenvalue 
equation for a discrete number of frequencies and then sorting by continuity of k. In contrast 
to the u(k) method, the periodicity is not relevant in this formulation of the problem and 
k(u) does not follow the first Brillouin zone. 

Because of the periodicity of the system, Bloch waves can be expanded in a series of har- 
monics where each harmonic corresponds with a value of k, if k is then a complex number, the 
evanescent behavior of a wave with a predetermined frequency would be multiexponential.— 
The complex band structures show the values of all of the complex values of k which con- 
tribute to the multi-exponential decay of the mode in the BG. As we will see later, for the 
case of the SC analyzed in this paper, we can only approximate the evanescent behavior in 
the modes inside the BG by considering the first term of this harmonic expansion in terms 
of k. 
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A. Supercell approximation 



One particularly interesting aspect of SCs is the possibility of creating point defects that 
confine acoustic waves in localized modes.— ^ Because of the locally breaking periodicity of 
the structure, defect modes can be created within the BG. These defect modes are strongly 
localized around the point defect: once the wave is inside the defect, it is trapped because the 
borders of the defect act as perfect mirrors for waves with frequencies in the BG. Localization 
depends on several parameters, such as the size of the point defect. However, in finite 
periodic structures the strength of sound localization also depends on the size of the structure 
around the defect because of the exponential decay of the outgoing wave.— 

To analyze the propagation of waves inside periodic structures with defects, authors 
have traditionally used PWE with supercell approximation. The supercell method requires 
the lowest possible interaction between defects. This results in a periodic arrangement of 
supercells that contain the point defect. With this method we can obtain the relation u(k) 
for crystals with local defects and, for instance, the physics of wave guides^ 8 or filters^ 
can be explained. 

In this section, we apply the supercell approximation to the EPWE. This methodology 
enables us to obtain the relation k(u) for defect modes. It will be interesting to discover 
how the imaginary part of the wave vector inside the BG changes with the creation of the 
defect. 

Consider an SC with primitive lattice vectors <5j (i = 1,2,3). The supercell is a cluster 
of rii x ri2 x 71% scatterers periodically placed in space. The primitive lattice vectors in the 
supercell approximation are a\ = n^i, and the complete set of lattices in the supercell 
approximation is {R'\R' = ka'i}, where and U are integers. The primitive reciprocal 
vectors are then 

g. = 2tt ^ ljka J X (20) 
a'l ■ (a' 2 x a' 3 ) 

where is the completely anti-symmetrical three-dimensional Levi-Civita symbol. The 
complete set of reciprocal lattice vectors in the supercell is {G\Gi = Nib'i} where N{ are 
integers. 

Finally, the structural factor of the supercell in this approximation has to be computed 
while taking into account the size of the supercell. If we consider a 2D SC with cylindrical 
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scatterers with a radius r and an n\ x ri2 sized supercell, the structure factor of the supercell 
is expressed by 

(ni-l)/2 (n 2 -l)/2 

jr(Q} = ^ ^ e < ia \eA+icL\G2\)P(G) ^21) 

i=-(n 1 -l)/2j=-(n 2 -l)/2 

where 

P(G) = %,UG). (22) 

/ is the filling fraction of the supercell, G — \G\ and a is the lattice constant of the 2D 
periodic system. 

By introducing the previous expressions in the matrices of the PWE ([8]), or in the case of 
the EPWE (fill , we can then use the supercell approximation to calculate the band structure 
of a periodic structure with, and without, a point defect. 



IV. NUMERICAL RESULTS 



We consider a 2D SC consisting of PVC cylinders of radius r in an air background arranged 
in a square lattice with a lattice constant a. The material parameters employed in the 
calculations are p air = 1.23kg/m 3 , ppvc = 1400kg/m 3 , c a i r = 340m/s and cpvc = 2380m/s. 
We consider a filling fraction / = 7tr 2 /a 2 ~ 0.65. We have used reduced magnitudes,- so 
the reduced frequency is Q = wa/(27rc/ losf ), and the reduced wave vector is K = ka/ (2tt). 



A. Complete array 

In Figure [T] we can observe the complex band structure obtained by EPWE for the SC 
described above. In the left panel we have represented the imaginary part of the wave 
vector in the TX direction; in the right panel we have shown the complex band structures 
in the TM direction; and the central panel shows the real part of the band structures. The 
imaginary part is not restricted in values of k; while the real part is restricted to the first 
Brillouin zone. The area in gray represents the full BG ranged between the frequencies 
Q\ = Uia/(27rChost) — 0.4057 and ^2 = W2fl/(27rc/ los j) = 0.7189. Note that the real part of 
the complex band structures has exactly the same values as in the case of the PWE. 

In Figure [1] we can observe that modes inside the BG present purely imaginary wave 
vectors and these can be characterized as evanescent modes with an exponential-like decay. 
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FIG. 1. (Color online) Band structure of an SC of PVC cylinders embedded in air with filling 
fraction / ~ 0.65. The left panel represents the imaginary part of the wave vector for each TX 
direction frequency. The central panel represents the real part of the wave vector, constrained in 
the first Brillouin zone, for each frequency. The right panel represents the imaginary part of the 
wave vector for each TM direction frequency. The red dashed line represents the imaginary part 
of the wave vector of the evanescent modes inside the BG. Reduced magnitudes have been used. 




The elegant and intuitive explanation of the evanescent behavior of modes inside the BG 
given by Joannopoulus^ is reproduced in Figured] in TX; as well as in TM directions (red 
dashed lines). The imaginary part of the wave number for frequencies inside the BG grows 
with values of frequency closer to the center of the BG; and disappears at the edges of the 
BG. In other words, the rate of decay is greater for frequencies closer to the center of the 
BG. We can also observe that the imaginary part of the wave vector connects propagating 
bands and so conserves the overall number of modes. 

A recent paper has shown the multi-exponential decay of evanescent modes in a photonic 
crystal.— In Figure El we can observe clearly that each frequency inside the BG is charac- 
terized by several values of Im(k), corresponding to the harmonics of the multi-exponential 
decay of the evanescent modes. In the Section [IV] we will see that only the first value of 
the Im(k) contributes to the decay of the mode, and therefore higher harmonics can be 
neglected and we can approximate in the same way as an exponential-like decay. 
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FIG. 2. (Color online) Band structure for an SC with an internal defect, calculated using the 
EPWE with supercell approximation. The left panel represents the imaginary part of the wave 
vector for each TX direction frequency. The right panel represents the real part, constrained in 
the first Brillouin zone, of the wave vector for each frequency. The green dashed line represents the 
frequency of the localized mode in the defect. The red dotted line represents the imaginary part 
of the wave vector of the evanescent modes inside the BG. Reduced magnitudes have been used. 

B. Defect modes 

In this paper, point defects have been created by removing cylinders in an SC. We 
have used the EPWE method with supercell approximation to analyze the propagating 
and evanescent behavior of modes in an SC with point defects. 

Figure |2] shows the complex band structures for the TX direction and real band structures 
for an SC with a point defect. In our case, we use only one direction of incidence to analyze 
the complex band structure because the localized mode appears at the same frequency for 
all the incidence directions. The supercell used for the calculations is shown in the inset of 
Figure [2j We can observe that the localized mode appears at Q 3 = u^a/ (2tt Ch os t) = 0.59 
(green dashed line). For frequencies in the BG, the borders of the point defect act as perfect 
mirrors and produce the localized mode in this cavity. The complex value of the k number 
for the modes inside the BG can be obtained by EPWE and becomes a purely real value 
for the localized mode (red dotted line and green dashed line). The value exactly coincides 
with the value obtained by PWE with supercell approximation. 
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V. EXPERIMENTAL RESULTS 



We performed the experiments in an echo- free chamber sized 8 x 6 x 3m 3 . To obtain the 
experimental dependence of the pressure all along the SC, we measured the pressure field at 
several points between two rows of the SC. To achieve this we built a finite SC and placed 
the microphone inside the periodic structure in a space between two rows. The finite 2D 
SC used in this paper was made of PVC cylinders hung in a frame and measuring 5ax5a. 
The radius of the cylinders was r = 10cm, and the lattice constant of the SC was a = 22cm. 
With these parameters, the finite SC has the same filling fraction (/ ~ 0.65) as in Section 
IIVI and the dimensions are large enough for the microphone to be placed between the rows. 
The microphone used was a prepolarized free-field 1/2" Type 4189 B&K. The diameter of 
the microphone was 1.32cm, which is approximately 0.06a, and so a low level of influence 
over the pressure field measured is expected. 

The 3DReAMS system is capable of sweeping the microphone through a 3D grid of 
measuring points located at any trajectory inside the echo-free chamber. The motion of 
the robot was controlled by an NI-PCI 7334. We analyzed the absolute value of the sound 
pressure between two rows of the SC by moving the microphone in steps of 1 cm. 

In Section [IV] we analyzed the upper and lower frequencies of the BG for an SC of PVC 
cylinders with the filling fraction value as in our experimental set up. By considering the 
corresponding values of the parameters of our experimental SC, we can obtain the frequency 
range of the BG. In our case, the BG appears between 627Hz and 1111Hz. To measure the 
propagation of sound inside the SC, we analyzed two different frequencies, one inside the 
BG and the other in the first transmission band. The frequencies were 920Hz and 442Hz, 
respectively. 

In Figure [3] we show the experimental measurements of the absolute value of the pressure 
inside SC for propagating and evanescent modes. These experimental results represents 
a novel measurement of the pressure field inside an SC. The inset of Figure [3] shows the 
measured points in steps of 1 cm placed between two rows of cylinders inside the SC using 
the 3DReAMS system. Blue squares with a continuous blue polygonal line represent the 
absolute value of the pressure of a frequency outside of the BG, that is 442Hz. This frequency 
represents a propagating mode inside the SC. Red circles with a polygonal red continuous 
line represent the absolute value of the pressure of a frequency inside the BG, that is 920Hz. 
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FIG. 3. (Color online) Absolute value of the pressure inside the SC in the positions between two 
rows. Blue squares represent these values for a frequency outside of the BG, 442Hz. Red circles 
represent these values for a frequency inside the BG, 920Hz. Black dots represent the values 
used to fit the exponential decay. Green line represents the fit of the exponential decay of the 
evanescent mode inside the BG. The black continuous line represents the absolute values of the 
pressure obtained by finite element methods. 

For the last case, we can observe the decay of the pressure inside the SC because of the 
evanescent behavior of the mode inside the BG. 

In contrast to the propagating mode (blue squares with a blue polygonal continuous 
line), the evanescent mode (red squares with a red polygonal continuous line) is practically 
extinguished at the end of the crystal - and just a small value remaining for the emerging 
pressure. This characteristic of evanescent behavior in finite SCs has been measured recently 
by Wu et al.— in an SC with a point defect. 

The value of the imaginary part of the first harmonic of the wave vector for the 920Hz 
frequency can be obtained from Figure HJ Using the values of parameters of the SC, we 
can observe a value Im(k) = — 5.6m _1 . From experimental data (see Figure |3]), we can 
fit the decay of the evanescent mode. We have chosen the points with maximum values in 
order to fit an exponential decay ae bx . The values of the fit are a = 0.05597 ± 0.0103Pa 
and b = Im(k) = —5.60 ± 1.45m _1 . Note that the experimental value is very close to the 
analytical value, i.e., the assumption that only the first harmonic is needed to represent the 
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multiexponential decay of the evanescent mode is correct. 

By solving the scattering problem inside the SC by means of the finite element method 
(FEM) we can analyze the evanescent behavior of the modes inside the BG of an SC. We 
have studied numerically the absolute value of the sound pressure between two rows of an 
SC. Continuity boundary conditions in the walls of the cylinders and the radiation condition 
at the borders of the numerical domain have been considered in the simulation. The black 
continuous line in Figure [3] represents the absolute values of pressure obtained numerically 
inside the SC, considering an incidence of a plane wave with a frequency of 920Hz. The 
correspondence between the experimental data (red polygonal line with open red circles) 
and the numerical results is clear. 

VI. CONCLUSIONS 

The propagation of waves inside periodic structures consists of propagating and evanes- 
cent modes. u(k) methods can be used to analyze the propagating modes, while evanes- 
cent modes are represented by the absence of k for some ranges of frequencies. In this 
paper, we extend the u(k) to the k(u>) method for the case of 2D SCs. We present the 
formulation of the supercell approximation for the k(u) method. For the EPWE we have 
predicted the evanescent nature of the modes inside the BG of an SC. In this paper we 
have reported measurements of the exponential-like decay of the acoustic field inside an SC. 
EPWE predicted a value for the imaginary part of the first harmonic of the wave number, 
Im(k) = —5.6m -1 ; and by fitting an exponential decay, ae bx , the experimental value we have 
obtained is b = Im(k) = —5.60 ± 1.45m _1 . Therefore, we can conclude that only the first 
harmonic contributes to the exponential-like decay of the evanescent mode. We have also 
shown that the imaginary part of the wave vector connects propagation bands and conserves 
the overall number of modes. 

We have also applied the EPWE with supercell aproximation to SC with point defects. 
We have analyzed the case of one vacancy observing the localized mode inside the BG 
predicted by EPWE. The value of the k number for this localized mode, that is purely 
imaginary in the case of complete SC, changes to purely real and it becomes in a passing 
mode as it was observed in the literature. The frequency of the localized mode exactly 
coincides with the value obtained by PWE. 
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Analytical, numerical, and experimental results reproduce with very good agreement the 
complex values of the wave vector inside the BG, meaning that these methodologies obtain 
good values for the exponential-like decay of the evanescent modes in an SC. This work 
shows the basis for the correct understanding of the design of narrow filters and wave guides 
based on phononic or sonic crystals with point defects. 
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